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ABSTRACT 

We use dynamical mass measurements of 16 black holes in transient low-mass X-ray binaries to infer 
the stellar black hole mass distribution in the parent population. We find that the observations are 
best described by a narrow mass distribution at 7.8 ±1.2 Mq. We identify a selection effect related 
to the choice of targets for optical follow-ups that results in a flux-limited sample. We demonstrate, 
however, that this selection effect does not introduce a bias in the observed distribution and cannot 
explain the absence of black holes in the 2 — 5 Mq mass range. On the high mass end, we argue that 
the rapid decline in the inferred distribution may be the result of the particular evolutionary channel 
followed by low-mass X-ray binaries. This is consistent with the presence of high-mass black holes 
in the persistent, high-mass X-ray binary sources. If the paucity of low-mass black holes is caused 
by a sudden decrease of the supernova explosion energy with increasing progenitor mass, this would 
have observable implications for ongoing transient surveys that target core-collapse supernovae. Our 
results also have significant implications for the calculation of event rates from the coalescence of black 
hole binaries for gravitational wave detectors. 
Subject headings: X-rays: binaries — black hole physics 



I. INTRODUCTION 

The distribution of stellar black hole masses in the Galaxy is intricately related to the population and evolution of 
massive stars, the energetics and dynamics of supernova explosions, and the dividing line between neutron stars and 
black holes. Inferring this distribution from observations helps address a number of outstanding questions in these 
areas. 

Observations of mass losing evolved stars are one of the direct ways of studying the end stages of stellar evolution 
(Smartt 2009). Similarly, supernova and transient surveys (e.g., Rau et al. 2009 and references therein) provide 
increasingly larger data samples that have revealed a range of supernova properties and energetics that is wider than 
previously anticipated. Studies of the end products of supernova explosions, namely of the black hole and neutron 
star mass distributions, provide a different and complementary approach to constraining models of massive stellar 
populations, evolution, and explosions (Timmes, Woosley, & Weaver 1996; Fryer 1999; Fryer et al. 2002; Woosley, 

; Heger, & Weaver 2002; Zhang, Woosley, & Heger 2008). 

The mass distribution of stellar black holes is a key ingredient in the calculation of the event rates for gravitational 
wave observatories such as LIGO, GEO500, and Virgo (Abadie et al. 2010). In addition, inferring this distribution 
may shed light on the origins and properties of pulsar-black hole binaries (Narayan, Piran, & Shemi 1991), which are 

' believed to be some of the best laboratories for strong- field gravity tests (Wex & Kopeikin 1999). 

During the last two decades, the masses of a large number of black holes in X-ray binaries have been measured 
(Remillard & McClintock 2006; McClintock & Remillard 2006), providing a sample suitable for the statistical inference 
of their parent mass distribution. Using this sample, we find in this article strong evidence for a very narrow distribution 
of masses for black holes in transient low-mass X-ray binaries. In particular, we show that there is a paucity of black 
holes with masses ~ 2 — 5 Mq. Our results confirm and strengthen an earlier finding of Bailyn et al. (1998), who 
argued for a low-mass gap based on a limited sample of black holes. We explore in detail below known selection biases 
and demonstrate that the low-mass gap cannot be attributed to the observational selection of targets. We further 
argue that the lack of high-mass black holes in low-mass X-ray binaries may result from a particular evolutionary path 
that leads to their formation. 

2. ACCRETING STELLAR-MASS BLACK HOLES AND BLACK HOLE CANDIDATES 

In Tables 1 and 2, we present basic data for all 23 confirmed black- hole X-ray binaries (Remillard & McClintock 
2006). For the purpose at hand, we divide these systems into three groups: eight transient systems with orbital periods 
exceeding one day; nine shorter-period transient systems for which the measurement of black hole mass is problematic; 
and six systems with persistent X-ray sources and massive 0/B-type secondaries. In these tables, we provide the 
following data for most of the systems: Galactic coordinates, maximum X-ray intensity, orbital period, a distance 
estimate, mass function, mass ratio, inclination angle, and black hole mass. Our primary focus is on the 17 transient 
systems listed in Tables 1 and 2. 

We summarize in Table 3 data for thirty-two additional transient systems that are believed to contain black hole 
primaries based on the spectral and timing properties of their X-ray sources. For each of these systems, we give both 
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celestial and Galactic coordinates and the maximum X-ray intensity that has been reported. 

2.1. The Twenty-Three Black Hole Binaries 

We now survey the constraints that have been placed on the masses of the 23 confirmed stellar black holes via 
dynamical measurements. In determining the masses of black holes in X-ray binaries, the mass function 

/(M) = = (1) 

is the most important and secure observable. The orbital period Porb and the half-amplitude of the velocity curve of the 
secondary star K can, in most cases, be determined precisely and accurately (Remillard & McClintock 2006; Charles 
& Coe 2006). These two quantities define the value of the mass function, which is an absolute lower limit on the mass 
of the compact object: i.e., M > ,f(M). The mass function relates the black hole mass M, the orbital inclination 
angle i, and the mass ratio q = M2/M, where M2 is the mass of the secondary star. Values of the mass function, the 
inclination, and the mass ratio for the 23 established black hole binaries are given in Table 2. The mass-ratio estimates 
were obtained or derived from Orosz (2003) and Narayan & McClintock (2005), with a few refinements based on the 
references cited in the table. Determining the black hole mass M, the quantity of interest, is challenging because in 
many cases it is difficult to obtain secure constraints on the inclination and the mass ratio. 

Table 1 separates the 23 systems into three groups based on their X-ray behavior: the persistent sources (PS), the 
long-period transients (LPT), and the short-period transients (SPT). This is important for understanding possible 
selection effects, which we explore in Section 4. Table 2, on the other hand, separates sources into three different 
groups based on the amount of data available on their mass ratios and inclinations. This impacts the inference of 
the masses of individual sources, which we discuss in Section 3. These tables sparely give for each source a few key 
references, which are supplemented by some additional references in the discussion that follows. 

2.1.1. The Persistent Sources 

These systems contain 0/B-type secondaries and are persistently X-ray bright. Five of them have relatively 
massive black-hole primaries, M > 10 Mq. However, the mass of LMC X-3 is presently poorly constrained: 
4 Mq < M < IIMq (Cowley 1992). In the case of Cyg X-1, a very wide range of masses down to 5 Mq 
has recently been considered (Caballero-Nievcs et al. 2009). However, recent VLB A observations have shown that 
the distance exceeds 1.5 kpc (M. Rcid, private communication), which firmly establishes M > 8 Mq. There are two 
caveats on the mass constraints on IC 10 X-1 and NGC 300-1: These results are based on less-reliable emission-line 
radial velocities and assume specific lower bounds on the masses of the secondary stars. 

2.1.2. The Long Period Transient Sources 

Nearly all transient black hole sources, which are fed by Roche-lobe overflow, exhibit long periods of deep quiescence 

during which the spectra of their secondary stars are prominent. For all eight of these long-period systems (Porb > 1 d) 
there is strong evidence that the masses of their black hole primaries exceed 6 Mq . We comment on the two weakest 
cases. (1) GX 339-4 has never reached a deep enough quiescent state to reveal its photospheric absorption lines, and 
its mass function was determined via the Bowen emission lines (Hynes et al. 2003). Furthermore, this system does not 
exhibit ellipsoidal light curves that allow its inclination to be constrained. For a defense of the M = 6 Mq lower limit 
see Munoz-Darias et al. (2008). (2) GRS 1915-1-105: The pioneering mass measurement of this system by Greiner et 
al. (2001) requires conflrmation. Their spectroscopic orbital period disagrees by 8% with a more recent photometric 
determination of the period (Neil et al. 2007). Again, there are no well-behaved ellipsoidal light curves that can be 
used to constrain the inclination, which in this case is inferred from a kinematic model of the relativistic radio jets 
(Mirabel & Rodriguez 1999). 

2.1.3. The Short Period Transient Sources 

Significant constraints have been placed on the black hole primaries of only three of the nine short period systems 
(Porb < 1 d). It is difficult to obtain reliable inclination constraints for these systems because studies in quiescence of 
their small, late-type secondaries are compromised by the presence of a relatively strong and variable component of 
non-stellar light (e.g., Zurita et al. 2003; Cantrell et al. 2008), which is continually present in all of these systems. 

A measurement of the inclination angle has been obtained only for A0620-00, the prototype system, which is two 
magnitudes brighter than the other short period transients. The study by Cantrell et al. (2010), a tour de force 
that makes use of 32 photometric data sets spanning 30 years, relies on 10 data sets obtained when the source was 
in a "passive" quiescent state (Cantrell et al. 2008). As indicated in Table 2, the authors constrain the mass to 
be M = 6.6 ± 0.25 Mq. Meanwhile, the masses of three other systems (XTE J1118-h480, Nova Oph 1977 and GS 
2000+251) have by dint of their large mass functions been shown to exceed 6 Mq (Table 2). The large mass function 
of XTE J1859+226 (Table 2) suggests that it also should be inchided in this group; however, this result is imreliable, 
having only been presented in an lAU Circular (Filippenko & Silverman 2001), and the orbital period is uncertain 
(Zurita et al. 2002). We exclude this source from further considerations on deriving the black-hole mass distribution. 

Among the 23 systems in Table 2, arguably the best candidate for hosting a low-mass black hole is GRO J0422-I-32. 
However, by briefly discussing recent attempts to determine the system's inclination and mass, we show that current 
results are presently unreliable: Beekman et al. (1997) constrained the inclination to lie in the range i = 10° — 31° 
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TABLE 1 

Properties of Twenty-Three Black Hole Binaries 





Coordinate 


Common Name 


Type'' 


I 


b 


Max. Int. 


^orb 


D 


References'^ 




Name 


or Prefix" 


(deg) 


(deg) 


(Crab) 


(hr) 


(kpc) 


(distance) 


1 


1354-64 


(GS) 


LPT 


310.0 


-2.8 


0.12 


61.1 


> 25 


1 


2 


1543-47 


(4U) 


LPT 


330.9 


+5.4 


15 


26.8 


7.5 ± 0.5 


2 


3 


1550-564 


(XTEJ) 


LPT 


325.9 


-1.8 


7.0 


37.0 


4.4 ± 0.5 


3 


4 


1655-40 


(GROJ) 


LPT 


345.0 


+2.5 


3.9 


62.9 


3.2 ± 0.5 


4 


5 


1659-487 


"V^ O O O A 


LPT 


338.9 


-4.3 


1.1 


42.1 


9 ± 3 


5 


D 






T PT^ 


D.o 




1 

±o 


O i .D 


y.y m z.4 


D 


7 


1 n 1 1 1 n 


(LtKo ) 


T OT^ 

l^r i 


45.4 


-0.2 


3. / 


739 


9 + 3 


'7 

7 


8 


o n o o 1 o o o 

2023+338 


(GS) 


LPT 


73.1 


-2.1 


20 


155.3 


2.39 ± U.14 


8 


9 


0422-1-32 


(GROJ) 


SPT 


166.0 


-12.0 


3 


5.1 


2 + 1 


9,10 


10 


0620-003 


(A) 


SPT 


210.0 


-6.5 


50 


7.8 


1.06 ± 0.12 


11 


11 


1009-45 


(GRS) 


SPT 


275.9 


+9.4 


0.8 


6.8 


3.82 ± 0.27 


10 






(^A i rjj ) 


or i 


10 t .0 


+ 02. 


U.U4 


A 1 

4. i 


1 v _L n 1 
L.i it U.i 


12 


13 


1124-683 


Nova Mus 91 


SPT 


295.3 


-7.1 


3 


10.4 


5.89 ± 0.26 


10 


14 


1650-500 


(XTEJ) 


SPT 


336.7 


-3.4 


0.6 


7.7 


2.6 + 0.7 


13 


15 


1705-250 


Nova Oph 77 


SPT 


358.2 


+9.1 


3.6 


12.5 


8.6 + 2.1 


14 


16 


1859-f226 


(XTEJ) 


SPT 


54.1 


+8.6 


1.5 


9.2^ 


8 + 3 


10 


17 


2000+251 


(GS) 


SPT 


63.4 


-3.0 


11 


8.3 


2.7 + 0.7 


14 


18 


0020+593 


IC 10 X-1 


PS 






0.00006 


34.9 






19 


0055-377 


NGC 300-1 


PS 






0.00004 


32.3 






20 


0133+305 


M33 X-7 


PS 






0.00002 


82.9 






21 


0538-641 


LMC X-3 


PS 






0.06 


40.9 






22 


0540-697 


LMC X-1 


PS 






0.03 


93.8 






23 


1956+350 


Cyg X-1 


PS 






2.3 


134.4 







" The entries in parentheses are prefixes to the coordinate names that identify the discovery X-ray mission. 
^ PS = persistent source; LPT = long period transient; SPT = short period transient. 

^ References: 1. Casares et al. 2009; 2. J. Orosz, private communication; 3. Orosz et al. 2010; 4. Hjellming & Rupen 1995; 5. 
Hynes et al. 2004; 6. Orosz et al. 2001; 7. Fender et al. 1999; 8. Miller-Jones et al. 2009; 9. Webb et al. 2000; 10. Hynes et al. 

2005; 11. CantrcU et al. 2010; 12. Gclino ct al. 2006; 13. Homan et al. 2006; 14. Barret et al. 1996 
^ Unconfirmed and uncertain value: sec text. 

and concluded M > 9 Mq. Webb et al. (2000) concluded i < 45° and M > 2.2 Mq with a maximum mass that is 
not constrained. Gelino and Harrison (2003) present the strongest evidence for a low mass; they conclude i > 43° and 
M < 4.92 Mq. However, in a more recent Keck ii'-band study Reynolds et al. (2007) found that "No clear ellipsoidal 
modulation is present in the light curve..." and concluded, "...that previous infrared-based attempts to constrain the 
mass of the putative black hole in this system are prone to considerable uncertainty." Thus, it appears that a far more 
comprehensive photometric study (cf. Cantrell et al. 2010) is required in order to obtain a firm mass constraint. 

There are six systems in Table 2 with black holes of indeterminate mass. The prospects for measuring or usefully 
constraining the masses of two of them, LMC X-3 and XTE J1859-226, are bright: Work on the former is almost 
complete (J. Orosz, private communication), and the additional data required to confirm the orbital period and large 
mass function of the latter are obtainable. That leaves four systems (GRO J0422+32, GRS 1009-45, Nova Mus 1991 
and XTE J1650-500) with black hole masses that are only very weakly constrained: M > f{M) > 1 — 3 Mq. As 
indicated above and as Cantrell et al. (2010) have shown, it will be very challenging to place stronger and reliable 
constraints on the masses of these black holes. 



2.2. Thirty-two Transient Black-Hole Candidates 

Data for 32 X-ray transient systems are given in Table 3. These systems lack radial velocity data, and most even lack 
an optical counterpart. Thus, presently, there are no dynamical constraints on the masses of their compact primaries, 
which are believed to be black holes because they share certain characteristic X-ray properties with the 23 established 
black holes (McClintock & Rcmillard 2006). As indicated in Table 3, the primary source of information about these 
systems is the catalogue of Liu ct al. (2007) and references therein. For additional information and references on many 
of these systems, see Table 4.3 and text in McClintock and Remillard (2006). 

3. BLACK HOLE MASS MEASUREMENTS AND CONSTRAINTS 

In this section, we use the measurements of the mass functions, as well as any available constraints on the mass 
ratios and inclinations for the black holes in low-mass X-ray binaries shown in Table 2 in order to place quantitative 
constraints on the individual black hole masses. In particular, our aim is to derive the likelihood Pj(data|M), which 
measure the chance of obtaining the particular set of data shown in Table 2 for the i-th source if that source had mass 
M. 

We divide the sources into three categories based on the amount and quality of information regarding their mass 
ratios and inclinations: 

(i) For six sources, the mass ratios and the inclinations are tightly constrained, leading to well-determined black hole 
masses. In this case, the probability distribution can be described as a Gaussian 



Pi(data|M) = C^exp 



(2) 
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with a mean Mo^i and a standard deviation aM,i- In this and the following expressions, Ci is a proper normalization 
constant such that 

/>oo 

Pi{d&ta\M)dM = 1. (3) 





This category includes A0620-003, 4U 1543-47, XTE J1550-564, GRO J1655-40, V4641 Sgr, and GS 2023+338. 
(ii) For the sources in the second category, there is only a measurement of the mass function and constraints on the 
mass ratio q. Here, we assume a Gaussian probability distribution over the mass function with a mean /o,i and a 
standard deviation cro,i. For the mass ratio, we adopt a uniform distribution 

P{q)dq = ^ (4) 

qmax Q'min 

between the minimum and maximum allowed mass ratios, q„iin and (/max, respectively. For each value of the mass 
ratio, the lack of eclipses implies a maximum value of the inclination, i.e., a minimum value of cosi, such that 

/ \ 1/3 

(cosi)„,in = 0.462 f^^j . (5) 
Assuming a uniform distribution over cosi subject to this constraint, i.e., 

I \ ,/ N d(cOSi) / N 1 /r.\ 

F{cosi\q)a{cosi) = r , (COSljmin < cosz < 1, (6) 

1 - (cOSZjmin 



yields 



Amin i(cosi)„„„ 1 - (cOSl)min [ 2aj^^ J 

The following nine sources belong to this category: GROJ 0422+32, GRS 1009-45, XTE J1118+480, Nova Mus 91, 
MS 1354-64, XTE J1650-500, GX 339-4, Nova Oph 77, and GS 2000+251. 

(Hi) This last category includes only GRS 1915+105, for which the mass function and the inclination have been 
measured, and the mass ratio has been constrained (see the discussion in Section 2 about the inclination measurement). 
In this case, we calculate the probability distribution over mass using equation ([T]), supplemented by a Gaussian 
distribution over inclination 



PM + un rl A I '^(cos») j [/o,,-Msin3^/(l + g)^]^ {i - ipf y 
P.(datalAf) = / dq 1^^^^^^^^ i„(eosz)..„ \ 2^^ 2^ ^ ' 



Figure [T] shows the likelihoods -Pi(data|M) for the 16 sources in the above three categories. The top panel includes 
sources in categories (i) and (Hi), while the bottom panel shows those in category (ii). 

A clustering of the observed black hole masses between ~ 6 — 10 Mq is already evident from Figure [TJ In the next 
section, we will carry out a formal Bayesian analysis to determine the parameters of the underlying mass distribution 
that is consistent with the observed Pi (data |M) shown here. If the likelihood for each source was narrow enough 
such that there was little or no overlap between them, then adding the likelihoods for the entire sample and coarsely 
binning the resulting distribution would provide a good estimate of the underlying mass distribution. Even though 
this condition is not entirely satisfied here, especially at the high mass end, we nevertheless show in Figure [2] this 
approximate mass distribution to get a sense of its gross properties. 

4. THE INTRINSIC DISTRIBUTION OF BLACK HOLE MASSES 

In this section, we use a parametric form of the black-hole mass distribution and the data discussed in Sections 2 
and 3 in order to determine its parameters. We will first consider an exponentially decaying mass distribution with a 
cut-off given by 

P(M-M ^ _ exp(Mc/M,cale) / exp(-M/M,eaic) , M>M, 

l-[M,M,,,i,,M,) - ^^^^^^ M<M, ■ 

This choice of the mass distribution is motivated by theoretical expectations based on the energetics of supernova 
explosions, as well as the density profiles and mass distributions of pre-supernova stars. The typical value of the mass 
scale is expected to lie in the range Mgcaie ~ 5.5 — 9 Mq (as we infer from the various figures in Fryer & Kalogera 
2001), whereas the cutoff mass is simply expected to be the maximum neutron star mass. Our goal is to find the 
values of the mass scale Mgcaio in the exponential and the cut-off mass that maximize a properly defined likelihood 
and to estimate their uncertainties. We will show below that the particular choice of the functional form of the mass 
distribution does not affect the main conclusions of the paper. 

In Section 3, we calculated, for each observed black hole, the probability Pi(data|Af), which measures the chance 
of making a particular observation if the black hole has mass M. What we want to calculate here is the probability 
P{Mscede, Mc\da.ta), which measures the likelihood of the parameters of the black hole mass distribution, given the 
observations. Using Bayes' theorem, we can write this as 

P(A/,eaic, Afcldata) - C2P(data|M,eaic, M,)P{M,,,i,)P{AQ , (10) 
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Fig. 1. — The likelihoods Pi(data|M) for the 16 sources in low-mass X-ray binaries that have been securely identified as black holes. The 
top panel includes sources in categories (i) and (in), while the bottom panel shows those in category (ii). The categories are based on the 
amount of information available on the mass ratios and inclinations of the black hole binaries and are discussed in more detail in the text. 

where C2 is the normahzation constant and P(Mscaic) and P(Mc) are the priors over the values of the mass scale and 
cut-off mass. We assume a flat prior over the mass scale between Mgcaie = and a maximum value Mgcaie — -^max, 
i.e., 

' 0, M,ealc < 

P(Mscalc) = { M^, < M,calo < M„,ax (11) 



0, 



The upper limit Mmax is imposed mostly for computational reasons and does not affect the results. We also adopt a 
similar prior over the cut-off mass between the maximum neutron-star mass, which we set to 2 Mq, and the minimum 
well-established mass measurement for a black hole. As will be evident from the results, the particular choice of this 
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Fig. 2. — The solid line shows the sum of likelihoods for the mass measurements of the 16 black holes in low- mass X-ray binaries. Note 
that because of the high-mass wings of the individual likelihoods, the shape of their sum is artificial at the high mass end. The dashed and 
dotted lines show the exponential and Gaussian distributions, respectively, with parameters that best fit the data (see §4). 



TABLE 2 

Dynamical Data for Twenty-Three Black Hole Binaries 





Coordinate 


Common Name 


fiM) 


Q 


i 


M 


References'^ 




Name 


or Prefix'' 


{Mq) 


^miii (/max 


(deg) 


(Mq) 




1 


0620-003 


(A) 


2.76 + 0.01 


0.056-0.064 


51.0 + 0.9 


6.6 + 0.25 


1,2 


2 


1543-47 


(4U) 


0.25 + 0.01 


0.25-0.31 


20.7+1.5 


9.4 + 1.0 


3 


3 


1550-564 


(XTEJ) 


7.73 + 0.40 


0.0-0.040 


74.7 + 3.8 


9.1 + 0.6 


4 


4 


1655-40 


(GROJ) 


2.73 + 0.09 


0.37-0.42 


70.2 + 1.9 


6.3 + 0.27 


3,5 


5 


1819.3-2525 


V4641 Sgr 


3.13 + 0.13 


0.42-0.45 


75 + 2 


7.1 + 0.3 


3 


6 


2023+338 


(GS) 


6.08 + 0.06 


0.056-0.063 


55 + 4 


12 + 2 


6 


7 


0422+32 


(GROJ) 


1.19 + 0.02 


0.076-0.31 






3,7 


8 


1009-45 


(GRS) 


3.17 + 0.12 


0.12-0.16 






3,8 


9 


1118+480 


(XTEJ) 


6.1 + 0.3 


0.035-0.044 






3,9,10 


10 


1124-683 


Nova Mus 91 


3.01 + 0.15 


0.11-0.21 






3,11 


11 


1354-64 


(GS) 


5.73 + 0.29 


0.08-0.15 






12 


12 


1650-500 


(XTEJ) 


2.73 + 0.56 


0.0-0.2'= 






13 


13 


1659-487 


GX 339-4 


5.8 + 0.5 


0.0-0.4'= 






14,15 


14 


1705-250 


Nova Oph 77 


4.86 + 0.13 


0.0-0.053 






3,6 


15 


2000+251 


(GS) 


5.01 + 0.12 


0.035-0.053 






3,6 


16 


1915+105 


(GRS) 


9.5 + 3.0 


0.025-0.091 


66+2 




16,17,18,19 


17 


1859+226 


(XTEJ) 


7.4+ 1.1^ 








20,21 


18 


0020+593 


IC 10 X-1 


7.64 + 1.26 






> 20 


22,23 


19 


0055-377 


NGC 300-1 


2.6 + 0.3 






> 10 


24 


20 


0133+305 


M33 X-7 


0.46 + 0.07 






15.65 + 1.45 


25 


21 


0538-641 


LMC X-3 


2.3 + 0.3 








26 


22 


0540-697 


LMC X-1 


0.886 + 0.037 






10.91 + 1.54 


27 


23 


1956+350 


Cyg X-1 


0.251 + 0.007 






> 8 


28,29 



^ The entries in parentheses arc prefixes to the coordinate names that identify the discovery X-ray mission. 
^ References: 1. Ncilscn ct al. 2008; 2. Cantrell et al. 2010; 3. Orosz 2003; 4. Orosz et aL 2010; 5. Greene et al. 2001; 6. 
Charles & Coe 2006; 7. Filipponko ct al. 1995; 8. Filippenko et al. 1999; 9. McClintock et al. 2001; 10. Gelino et al. 2006; 
11. Orosz et al. 1996; 12. Casaros et al. 2009; 13. Orosz et al. 2004; 14. Hynes et al. 2003; 15. Muiioz-Darias et al. 2005; 16. 
Greiner et al. 2001; 17. Neil ct al. 2007; 18. Harlaftis & Greiner 2004; 19. Fender et al. 1999; 20. Filippenko & Chornock 
2001; 21. Zuritaet al. 2002; 22. Prestwich et al. 2007; 23. Silverman & Filippenko 2008; 24. Crowther et al. 2010; 25. Orosz 
ct al. 2007; 26. Cowley 1992; 27. Orosz et al. 2009; 28. Caballero-Nieves et al. 2009; 29. M. Reid, private communication. 
^ Estimated range based on extreme values observed for systems with comparable orbital periods. 
^ Unconfirmed and uncertain value; see text. 
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Fig. 3. — The parameters of an exponential black hole mass distribution with a low-mass cutoff. The cut-off mass is well above theoretical 
expectations, indicating a sizable gap between neutron-star and black hole masses. Furthermore, the mass scale in the exponential is 
significantly smaller than theoretical expectations. 

range does not affect the measured parameters. We also repeated this analysis with logarithmic priors in the two 
parameters and found that the results are insensistive to the choice of priors. 

In equation (|10p . the quantity P(data|Mscaiej Mc) measures the chance that we make a particular set of observations 
for the ensemble of black holes, given the values of the parameters of the mass distribution. We need now to estimate 
this quantity, given the likelihoods for the individual sources. We will assume that each measurement is independent 
of all the others, so that 

P(data|M,caic, A/c) = 11/ dMP,{d&ta\M)P{M ; M.eaic, M,) . (12) 
Combining this last equation with equation (jlOp we obtain 

P(Afscaic, Mcldata) = CP(Af,caic)P(Mc) H / dM P,{data\M)P{M; M.^aie, M,) , (13) 

i 

where C is the overall normalization constant. 

We show in Figure [3] the 68% and 95% confidence contours of the mass scale and cut-off mass that best describe 
the observations and compare them to the theoretical expectation. The lack of black holes below ^ 5 Mq and the 
rapid decline of the exponential distribution at the high mass end are both remarkable (see the dashed line in Fig. [5]). 
The latter result is not at odds with the relatively high mass of GRS 1915-1-105 because of the wide and shallow mass 
probability distribution of this source. 

In order to explore whether the small number of sources with very well-determined masses dominate this result, we 
did the following test. We repeated the calculation using the mass functions and constraints on the mass ratios for all 
sources but ignoring any information on the inclinations of the binary systems. This, in effect, is equivalent to treating 
all sources in category (i) using the formalism we applied to sources in category (ii), integrating over all possible values 
of inclination. Figure 0] shows the 68% and 95% confidence contours of the parameters of the exponential distribution 
in this test case. Although the allowed range of values is increased, as expected, the low-mass gap and the discrepancy 
with the theoretically expected mass scale remain robust. 

The narrowness of the mass distribution implied by the above results motivated us to explore different functional 
forms of the underlying distribution, and in particular, a Gaussian function. The Gaussian function here serves as a 
phcnomcnological two-parameter description of a narrow distribution and is not necessarily motivated by theory. We 
show in Figure [5] the parameters of such a Gaussian distribution that best describes the observations. The masses 
of all 16 black holes are consistent with a narrow distribution at 7.8 ± 1.2 Mq. This result is in agreement with an 
earlier, more limited, study by Bailyn et al. (1998). 
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Fig. 4. — The parameters of an exponential black- hole mass distribution when all possible inclinations for each binary system were 
allowed (i.e., no inclination information was taken into account). The existence of the low-mass gap and the small exponential mass scale 
remain robust. 




Fig. 5.— 



The parameters of a Gaussian black hole mass distribution. The data are best described by a narrow distribution at 7.8±1.2 Mq. 



9 



TABLE 3 

Thirty-Two Transient Black-Hole Candidates 



Source Name RA{J2000) Dec{J2000) i(deg) 6(deg) I(Crab)'' References'" 



I 


1 A 1 524-61 


15 


28 


17.2 


-61 


52 


58 


320.3 


-4.4 


0.95 


I 


2 


4TI 1630-47 


16 


34 


01.61 


-47 


23 


34.8 


336.9 


+0.3 


1.4 


1 


3 


XTE J1652-453 


16 


52 


20.33 


-45 


20 


39.6 


340.5 


-0.8 


0.19 


2 


4 


TOR 117098-3628 


17 


09 


45.93 


-36 


27 


58.2 


349.6 


+2.1 


0.13 


1 


5 


SAX 11 71 1 6-3808 


17 


11 


37.1 


-38 


07 


05.7 


348.6 


+0.8 


0.13 


I 


g 


GRO 11719-24 


17 


19 


36.93 


-25 


01 


03.4 


359.9 


+7.0 


1.5 


1 


7 


XTE J1720-318 


17 


19 


58.994 


-31 


45 


01.25 


354.6 


+3.1 


0.41 


1 


g 


lOR 117269-4737 


17 


26 


49.28 


-47 


38 


24.9 


342.2 


-6.9 


0.083 


I 


9 


HRS 1 730-31 2 


17 


33 


52.3 


-31 


12 


25 


356.7 


-H 


0.20 


I 


10 


HRS 1 737-31 


17 


40 


09 


-31 


02 


24 


357.3 


+0.6 


0.026 


I 


11 


HRS 1 739-278 


17 


42 


40.03 


-27 


44 


52.7 


0.7 


-H 2 


0.86 


I 


12 


KS 11739-304 


17 


42 


44.0 


-30 


30 


51 


358.3 


-0.3 


0.009 


1 


13 


1 F, 1 740 7-2942 


17 


43 


54.83 


-29 


44 


42.6 


359.1 


-0.1 


0.03 


I 


14 


1 A 1 742-289 


17 


45 


37.0 


-29 


01 


07 


359.9 


-0.0 


2.0 


I 


15 


HI 743-322 


17 


46 


15.57 


-32 


14 


01.1 


357.1 


-1.6 


0.77 


I 


16 


YTF, T1 74S 9S)R 

_L LLi O 1- 1 rlO ZlOO 


17 


48 


05.06 


-28 


28 


25.8 


0.7 


-0.2 


0.64 




17 


lOR 117497-2821 


17 


49 


38.037 


-28 


21 


17.37 


1.0 


-0.5 


0.03 


I 


18 


ST,X 1 746-331 


17 


49 


48.3 


-33 


12 


26 


356.9 


-3.1 


0.27 


I 


19 


Swift J1752-223 


17 


52 


15.10 


-22 


20 


32.78 


6.4 


+2.1 


0.11 


3 


20 


Swift J1753.5-0127 


17 


53 


28.29 


-01 


27 


06.22 


24.9 


+12.2 


0.16 


4 


21 


XTE J1755-324 


17 


55 


28.6 


-32 


28 


39 


358.0 


-3.6 


0.18 


1 


22 


4U 1755-33 


17 


58 


40.0 


-33 


48 


27 


357.2 


-4.9 


0.10 


1 


23 


GRS 1758-258 


18 


01 


12.40 


-25 


44 


36.1 


4.5 


-1.4 


0.02 


1 


24 


SAX J1805. 5-2031 


18 


05 


34 


-20 


30 


48 


9.6 


+0.3 


0.051 


1 


25 


XTE J1817-330 


18 


17 


43.54 


-33 


01 


07.8 


359.8 


-8.0 


1.35 


1 


26 


XTE J1818-245 


18 


18 


24.4 


-24 


32 


18.0 


7.4 


-4.2 


0.51 


1 


27 


Swift J1842.5-1124 


18 


42 


17.33 


-11 


25 


00.6 


21.7 


-3.2 


0.04 


5 


28 


EXO 1846-031 


18 


46 


39.8 


-03 


07 


21 


29.9 


-0.9 


0.3 


1 


29 


IGR J18539+0727 


18 


53 


54 


-1-07 27 


39.8 


+2.8 


0.02 


1 


30 


XTE J1856+053 


18 


56 


39 


+05 19 48 


38.3 


+1.3 


0.07 


1 


31 


XTE J1908+094 


19 


08 


53.077 


+09 23 04.90 


43.3 


+0.4 


0.10 


1 


32 


XTE J2012+381 


20 


12 


37.71 


+38 11 01.1 


75.4 


+2.2 


0.16 


1 



^ Approximate maximum X-ray intensity from the Liu et al. (2007) catalogue supplemented in a few cases 
by data from other sources. In some instances, the intensity is for a higher enery band than the standard 
~ 2 - 10 keV band. 

^ References: 1. Liu et al. 2007; 2. Markwardt et al. 2009; 3. Torres et al. 2009: 4. Fender et al. 2005; 5. Krimm 
et al. 2008. 



5. OBSERVATIONAL SELECTION EFFECTS 

In the previous section, we discussed two distinct results we obtained on the mass distribution of stellar black holes 
using the data on 16 sources in low-mass X-ray binaries: The cut-off mass at the low end is > 5 Mq, indicating a 
significant lack of black holes in the ~ 2 — 5 Mq range, and the mass scale in the exponential is ~ 1.2 Mq, indicating a 
rapid decline of the distribution at the high-mass end. In other words, the black hole mass distribution is narrow and 
can also be represented with a Gaussian at 7.8 ± 1.2 Mq. We now explore two questions about our results: (i) whether 
any observational selection effects could have biased the black-hole mass distribution we inferred for transient sources 
and (ii) whether the black holes in transient low- mass X-ray binaries are a representative sample of the population of 
stellar black holes in the Galaxy. 

5.1. Selection Biases in the Mass Distribution of Black Holes in Transient X-ray Binaries 

The masses of stellar black holes have generally been measured when the black holes (i) reside in binary systems, 
(ii) are transient X-ray sources, and (Hi) reach peak fluxes above 1 Crab (about 2.5 x 10~^ erg s~^ cm~^) during their 
outbursts. The first requirement naturally arises from the nature of the dynamical mass measurements. The second 
requirement ensures that optical observations of the companion star can be carried out during the quiescent phase so 
that the radial velocity measurements are not contaminated by the emission from the accretion disk'^. Finally, the 
third requirement is a common strategy of observing campaigns that aims to select among the transients those sources 
that are the closest and the least obscured. 

Black hole transients have been observed in the last four decades with a large collection of X-ray telescopes. We, 
therefore, tried to verify and quantify the third requirement for follow-up by compiling the highest flux level observed 
from each transient black hole or black hole candidate in outburst by any instrument (see Table 1). We show in 
Figure [6] the histogram of the highest flux levels observed during outburst from the 16 confirmed black holes (in red) 
and from the black hole candidates (in blue). The dichotomy between the fluxes of the two groups is striking and 
indeed shows that primarily the brightest sources have been followed up for dynamical mass measurements. 

We then explored whether the 1 Crab flux threshold during outburst could lead to a selection bias against low-mass 
black holes. Setting a flux limit for optical follow-ups restricts, in principle, the volume of the Galaxy that is sampled. 

^ The dynamical Bowen technique (Steeghs & Casares 2002) offers the possibihty of measuring the masses of black holes in persistent 
systems. 
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Fig. 6. — The histogram of maximum flux observed during the outbursts of black holes with dynamical mass determinations (red) and 
other black hole candidates (blue). A clear dichotomy between the two groups at 1 Crab is evident. 





Fig. 7. — The positions on the Galactic disk of the black holes with secure dynamical mass measurements. The bakground grayscale 
image is the electron distribution in the Galaxy according to the model of Cordes &: Lazio (2001). The plus sign represents the position of 
the Sun. 
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If, in addition, the low-mass black holes form an intrinsically fainter population than that of more massive black holes, 
this could indeed introduce a selection bias. 

Consider, for example, the case in which the maximum luminosity during the outburst of an X-ray transient is 
proportional to the mass of the black hole, so that the binary systems with low-mass black holes are intrinsically 
dimmer than their high mass counterparts. Because the volume of the disk component of the Galaxy scales as the 
distance squared, the volume sampled by sources selected for follow-up observations scales as the mass of the black 
hole. If the mass distribution of black holes was rather flat, then, on average, ^ 4 times more ^ 8 M© black holes 
would be followed up than ^ 2 M© black holes, leading to an artificial reduction in the number of low-mass black 
holes. 

If there were indeed a population of intrinsically fainter, low-mass black holes, for the 16 more massive ones that have 
been observed, 4 low mass black holes should have been found within the same flux limits according to the previous 
argument. In contrast, there is none in the current sample. This fact argues that the lack of low-mass black holes is 
real. 

Another way to determine whether the volume-limited sample introduces a mass bias is to compare the number and 
distribution in the Galaxy of secure black hole sources to those of the black hole candidates. We plot in Figure [7] the 
galactic distribution of the black holes with secure mass determinations. It is evident that these sources lie within 
~ 1/3 of the volume of the galactic disk. As a result, if the observed black holes are a representative sample of the 
true population, our catalogue of black-hole candidates should contain twice as many sources that are fainter simply 
because they lie further away in the Galaxy compared to the 16 sources that have been followed up. This is precisely 
the ratio of the number of black-hole candidate sources in Table 3 to that of secure black holes in Table 1: 32:16. 

The 16 secure black holes have a median intensity of 3.6 Crab and a median distance of ^4 kpc. The 32 black 
hole candidates, on the other hand, have a median intensity of 0.18 Crab, so we expect them to be on average 
•^3. 6/0. 18 = 4.5 times farther. (Note, however, that confirmed black holes have benefited from almost continuous 
monitoring that led to an accurate measurement of their the maximum intensities. The black hole candidates have 
been observed more sparsely, so that the maximum reported intensity is in most cases likely to be below the actual 
maximum intensity that the sources reached). This further supports the idea that the black hole candidates lie in the 
far two-thirds of the Galaxy. 

We, therefore, conclude that there is no room within the 32 black-hole candidates to hide any sizable population of 
sources that are fainter because they are intrinsically less massive. The fact that our sample of secure black-holes is 
flux limited does not imply that it is luminosity limited, and, therefore, does not introduce any selection bias against 
the detection of low-mass black holes. 

5.2. Are Black Holes in Transient Low-Mass X-ray Binaries a Representative Sample? 

The sample of 16 sources that we used to infer the mass distribution of stellar black holes has been drawn from 
transient low-mass X-ray binaries^. Here, we consider the possibility that low-mass black holes exist, but not in 
transient low-mass X-ray binaries. In turn, we will consider the possibility that low-mass black holes are primarily in 
persistent binaries or that they are not members of binary systems at all. 

(i) Are accreting low-mass (2 — 5 Mq) or high mass f> 9 Mq) black holes persistent X-ray sources? — Most known 
neutron-star low-mass X-ray binaries are persistent sources, whereas the majority of black hole binaries are transient 
systems. (The latter statement is primarily a consequence of the fact that the transient behavior has so far been 
required in establishing the black hole nature of a compact object in a majority of the sources). Is it possible that 
low-mass black holes behave more like neutron stars than their more massive counterparts? 

The most prominent explanation of the transient behavior of X-ray binaries is the irradiated disk instability model 
(van Paradijs 1996; King, Kolb, & Burderi 1996; Dubus et al. 1999). In this model, whether an X-ray binary is a 
transient depends only weakly on the nature of the primary but is mainly determined by the orbital period of the 
system, the mass of the primary, and the evolutionary state of the companion star (King et al. 1996). Thus, the 
expectation is that a 4 Mq black hole binary should behave like one containing a 6 Mq black hole if they have 
comparable orbital periods. Furthermore, our sample of 16 transient black hole sources includes both short-period 
and long-period systems, indicating that the transient behavior is seen over a range of orbital periods. We, therefore, 
find it unlikely, albeit still possible, that low-mass black holes would preferentially exist in persistent binaries if the 
distribution of their orbital periods is similar to that of the confirmed black holes. 

We, nevertheless, explore whether the known number of persistent X-ray sources in the Galaxy is large enough to 
harbor the low-mass back holes if, indeed, all of the latter are persistent systems. To this end, we analyzed the catalogue 
of low-mass X-ray binaries by Liu, van den Heuvel, & van Paradijs (2007). Starting with the total number of non- 
pulsing X-ray binaries in the Galaxy, we identified the following categories: (a) Confirmed black holes, which refer to 
those with dynamical mass measurements; (b) Black hole candidates, which are sources with no mass determinations, 
but possess spectra and radio properties that resemble those of the confirmed black holes; (c) Unidentified sources, 
which are those with significant number of observations but no distinguishing characteristics to reveal the nature of 
the primary; (d) Neutron stars, which are sources that exhibit thermonuclear X-ray bursts; and (e) Sources for which 
there is very little information. In Table 4, we show the numbers and relative fractions of sources in these categories. 

Table 4 indicates that there are a total of 72 known persistent X-ray sources. 46 out of these 72 sources have shown 
thermonuclear bursts and can securely be identified as neutron stars, leaving at most 26 sources as possible low-mass 

^ Persistent sources with O/B-type secondaries have also been successfully targeted for mass measurements; see Tables 1 and 2. 
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black hole binaries. In reality, only two of these remaining 26 sources show evidence for black hole nature, and a 
further 7 show no thermonuclear bursts even though they have been extensively observed. We will, nevertheless, allow 
the possibility that all can harbor low-mass black holes. 

For any integrable (i.e., declining faster than logarithmic) distribution of black hole masses that does not have a 
low-mass gap, there should be more black holes with masses in the 2 — 5Mq range than with masses > 5Mq. If all 
the 26 persistent sourcse were low-mass black holes, we would expect this number to be larger than the total number 
of transient (higher mass) black hole sources in the Galaxy. From the galactic distribution of transient black-hole 
binaries, we estimated (see previous subsection) that there are at least three times as many black holes in the Galaxy 
with masses > 5Mq than our current sample of 16. Taking into account the short duty cycles of transient black- hole 
binaries, we can safely infer that the total number of such systems in the galaxy is » 48, which is significantly larger 
than 26. As a result, even if all the 26 non-bursting sources were low-mass black holes, they would not be enough to 
efficiently close the low-mass gap between ^ 2 and ^ 5 Mq in the inferred distribution. Nevertheless, applying the 
dynamical Bowen technique to persistent low-mass X-ray binaries will help explore further this possibility. 

(ii) Are low-mass (2 — 5 Mq) or high niMss f> 9Mq) black holes not members of X-ray binary systems? — The fraction 
of binaries that survive the supernova explosion is determined by a number of factors, such as the orbital separation, 
the amount of mass loss and the kick velocity during the supernova explosion. These factors depend on the progenitor 
and the remnant masses and may work in a way to hinder the formation of low-mass X-ray binaries with low- or 
high-mass black holes. Moreover, even if the binaries survive the explosions, systems with low- or high-mass remnants 
may become wider and perhaps never reach contact. It appears unlikely that this is the explanation of the paucity of 
low- mass black holes between the systems with dynamical mass measurements (sec Fryer & Kalogera 2001), although 
addressing this question would only be possible through detailed population synthesis calculations. 

On the other hand, the particular evolutionary path that leads to the formation of black- hole X-ray binaries might be 
responsible for the lack of systems with relatively massive black holes (> 9 Mq). Mass transfer in the low-mass X-ray 
binary phase requires small orbital separations between the black hole and the low-mass companion. The common 
envelope phase, which is invoked to reduce sufficiently the orbital separation of the binary before the supernova, is 
thought to lead to the expulsion of the hydrogen envelope of the pre-supcrnova star, leaving behind a bare helium core. 
Furthermore, the winds from the resulting helium cores in a Wolf-Rayet phase are expected to lead to further mass 
loss, albeit at rates that are highly uncertain. As a result, it is possible that the black hole masses in contact binaries 
are capped at ~ 10 Mq, although this value depends strongly on the assumed mass loss rates (Woosley, Heger, & 
Weaver 2002). This evolutionary path could provide a natural explanation for the rapid decline of the inferred mass 
function at the high-mass end. The same mechanism does not necessarily cap the mass of the black holes in high-mass 
X-ray binaries, which are wind-fed. 

6. IMPLICATIONS FOR BLACK HOLE FORMATION 

We now explore the implications of the narrow mass distribution of black holes in low-mass X-ray binaries for the 
progenitors of these systems and the supernovae that form the black holes. 

As we discussed in the previous section, the rapid decline at the high mass end of the inferred distribution is 
plausibly the result of the particular binary formation channel that leads to low-mass X-ray binaries. In contrast to 
this population, relatively more massive black holes (> 20 Mq) have been detected in high mass X-ray binaries (see 
the persistent sources in Table 2). This argues that there is no intrinsic barrier for the formation of massive black 
holes as the end stages of stellar evolution but that the differences between the black hole masses in the two binary 
populations are due to the different evolutionary paths that they follow. Observational studies of mass loss from stars 
with hydrogen and helium envelopes, in isolation and in tight binary systems, will elucidate and test these ideas. 

The low-mass gap, on the other hand, cannot be accounted for by observational selection effects (even though the 
mass measurement method requires a particular choice of targets and the observing strategy has resulted in a flux- 
limited sample). It is also unlikely that it results from the particular evolutionary paths of low-mass X-ray binaries. 
Therefore, it appears to be real. 

The observed distribution of masses could be different than the distribution of birth masses because mass transfer 
in the X-ray binary phase could have caused the black holes to evolve away from the masses they were born with. 
Population synthesis models (e.g., Fragos et al. 2010) show, however, that the least massive black holes are preferentially 
born with low-mass companions (< 2Mq) and, therefore, could not accrete a significant amount of matter over the 
lifetime of the binary. When the companions are more massive than 5 Mq at the start of the mass transfer phase, the 
black holes are born with masses that are also > 5 Mq. As a result, the accretion in the binary phase is not likely to 
bridge the gap between 2 and 5 Mq. 

It has been suggested that black hole evaporation in braneworld gravity models can lead to a deficiency of low-mass 
black holes in the steady state population. A gap could be created in this context because the rate of evaporation 
in braneworld gravity is very rapid and increases with decreasing black hole mass (Postnov & Cherepashchuk 2003). 
However, recent constraints on the rate of evaporation obtained using the current population of black holes preclude 
this possibility (Johannsen, Psaltis, & McClintock 2009). 

The paucity of black holes with masses less than 5 Mq are likely, therefore, to be related to the physics of supernova 
explosions that lead to the formation of black holes. Numerical simulations of supernova explosions typically generate 
a continuous distribution of black hole masses that decays as an exponential (Fryer 1999; Fryer & Kalogera 2001). 
The continuity of masses is primarily a consequence of the relatively gradual dependence of explosion energies on the 
masses of the progenitors and the fact that explosion energies are still comparable to (although smaller than) the 
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TABLE 4 

The Population of Low-Mass X-ray Binaries in 
THE Galaxy. 



Primary 


Type 


Number 


Fraction 


Neutron Star 


Persistent 


46 


28% 


Neutron Star 


Transient 


39 


23% 


Confirmed BH 


Persistent 





0% 


Confirmed BH 


Transient 


16 


9% 


BH Candidate 


Persistent 


2 


1% 


BH Candidate 


Transient 


30 


18% 


Unidentified 


Persistent 


7 


4% 


Unidentified 


Transient 


3 


2% 


Little Information 


Persistent 


17 


11% 


Little Information 


Transient 


7 


4% 



binding energy of the stellar envelopes. Although creating a mass gap is difBcult given the current understanding of 
the supernova energetics, it has been suggested that it can be achieved under the ad hoc assumption that the explosion 
energy has a step-function dependence on progenitor mass and that it plunges to zero for stars more massive than 
- 25 Mq (Fryer & Kalogera 2001). 

Such a bimodality in the energies of explosions that form neutron stars versus black holes should become apparent in 
the large sample of supernovae anticipated from the ongoing surveys that are sensitive to underluminous core-collapse 
supernovae. 
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